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ABSTRACT 

We present early results from an ongoing study of the kinematic structure of star-forming galaxies 
at redshift z ^ 2 ~ ?i using integral-field spectroscopy of rest-frame optical nebular emission lines 
in combination with Keck laser guide star adaptive optics (LGSAO). We show kinematic maps of 
3 target galaxies Q1623-BX453, Q0449-BX93, and DSF2237a-C2 located at redshifts z = 2.1820, 
2.0067, and 3.3172 respectively, each of which is well-resolved with a PSF measuring approximately 
0.11 - 0.15 arcsec (~ 900 - 1200 pc at z ~ 2 — 3) after cosmetic smoothing. Neither galaxy at z ~ 2 
exhibits substantial kinematic structure on scales > 30 km s~^; both are instead consistent with 
largely dispersion-dominated velocity fields with cr ~ 80 km s^^ along any given line of sight into the 
galaxy. While the primary emission component of Q0449-BX93 shows no spatially-resolved kinematic 
structure, a faint, kinematically distinct emission region is superposed on the primary region at a 
relative velocity of ^ 180 km s"^, suggesting the possible presence of a merging satellite galaxy. 
In contrast, DSF2237a-C2 presents a well-resolved gradient in velocity over a distance of ~ 4 kpc 
with peak-to-peak amplitude of 140 km s^^. This velocity shear was previously undetected in seeing- 
limited long-slit observations despite serendipitous alignment of the slit with the kinematic major axis, 
highlighting the importance of LGSAO for understanding velocity structure on subarcsecond scales. 
It is unlikely that DSF2237a-C2 represents a dynamically cold rotating disk of ionized gas as the 
local velocity dispersion of the galaxy {a = 79 km s~^) is comparable to the observed shear. Using 
extant multi-wavelength spectroscopy and photometry we relate these kinematic data to physical 
properties such as stellar mass, gas fraction, star formation rate, and outflow kinematics and consider 
the applicability of current galaxy formation models. While some gas cooling models reproduce 
the observed kinematics better than a simple rotating disk model, even these provide a poor overall 
description of the target galaxies, suggesting that our current understanding of gas cooling mechanisms 
in galaxies in the early universe is (at best) incomplete. 

Subject headings: galaxies: high-redshift - galaxies: kinematics and dynamics - galaxies: starburst 



1. INTRODUCTION 

At high redshifts z ~ 2 — 3 galaxies undergo strong 
and rapid evolution from highly irregular clumps of star- 
formation into the familiar Hubble sequence of the local 
universe (Giavalisco et al. 1996, Papovich et al. 2005), 
during which time they are thought to accumulate the 
majority of their stellar mass (Dickinson et al. 2003). 
The evolution of structure and ordered motion in galax- 
ies at such redshifts is therefore of considerable interest 
in the attempt to understand the process of galaxy for- 
mation. As such, numerous authors (e.g. Steidel et al. 
2004, Papovich et al. 2006, Reddy et al. 2006a, Erb et 
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al. 2006b, and references therein) have used deep UV-IR 
photometry and spectroscopy to ascertain such global 
properties as stellar mass, population age, AGN frac- 
tion, star formation rate, and the strength of large-scale 
gaseous outflows expelled by rapid starbursts. 

These studies have greatly increased our understand- 
ing of the physical nature of these galaxies, indicating 
(for instance) that the z ~ 2 universe contains galaxies 
in a wide variety of evolutionary states from young and 
actively star-forming to massive and passively evolving 
(Glazebrook et al. 2004, Daddi et al. 2005, Papovich 
et al. 2006, Reddy et al. 2006a, and references therein) 
and suggesting (based on clustering properties) that typ- 
ical star-forming galaxies at z ^ 2 may evolve into the 
elliptical galaxy population of the local universe (Adel- 
berger et al. 2005a). However, in general it is unknown 
whether the majority of star formation in these galaxies 
occurs in flattened disk-like systems (as anticipated by 
cosmological models, e.g. Mo, Mau, & White 1998) or 
disordered non-equilibrium configurations, and whether 
this star formation is due primarily to large-scale gravita- 
tional instabilities, tidal shocks induced by major merg- 
ers, or some other as-yet undetermined physical mecha- 
nism. Similarly, it is unknown whether individual regions 
of star formation within a given galaxy follow a global 
abundance pattern or exhibit strong variations in chem- 
ical enrichment. 
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These questions are most readily investigated by means 
of optical nebular emission lines (redshifted into the near- 
IR at z > 1) such as Ha, H/3, [Oiii] A5007, and [Nii] 
A6584 which predominantly trace ionized gas in active 
star forming regions and can provide a kinematic probe 
of the gravitational potential. Such studies using slit 
spectroscopy (e.g. Erb et al. 2003,2004,2006a; Weiner 
et al. 2006) have met with moderate success, enabling 
the determination of dynamical masses for a large sample 
of galaxies at z > 1, but are complicated by uncertain- 
ties resulting from slit placement, possible misalignment 
with the kinematic major axis, and the small angular 
size of the galaxies (typically less than an arcsecond; 
see Fig. 2 of Law et al. 2007) relative to the size of 
the seeing disk. The twin complications of slit place- 
ment and alignment have recently been overcome using 
seeing-limited integral-field unit (IFU) spectroscopy to 
obtain two-dimensional maps of the ionized gas kinemat- 
ics within these galaxies. A sample of the largest and 
brightest galaxies in the rest-frame UV-selected sample 
has recently been observed by Forster Schreiber et al. 
(2006), who find (in 9 out of 14 cases) smoothly vary- 
ing velocity shear generally consistent with rotation in a 
flattened disk-like configuration. However, it is unknown 
whether atmospheric seeing (typically ~ half the size of 
the galaxy) masks more complex structure, and to what 
extent the prevalence of velocity shear in these sources is 
a product of the selection criteria of the study, which 
(for observational reasons) preferentially select bright, 
extended sources many of which exhibited large veloc- 
ity shear or dispersion in previous long-slit studies (Erb 
et al. 2003, 2006a). 

With the aid of laser guide star adaptive optics 
(LGSAO) technology it is now possible to study high- 
redshift galaxies on sub-kiloparsec scales, permitting 
fruitful observation of galaxies either unresolved or 
barely resolved in previous ground-based observations. 
Wright et al. (2007) have demonstrated the utility of 
LGSAO in combination with integral-field spectroscopy 
for a galaxy at redshift z ~ 1.5; in this contribution we 
present early results from an ongoing study of galaxies in 
the redshift range z ~ 2 — 3. In §2 we describe our obser- 
vational program and basic data reduction techniques. 
In §3 and 4 we describe the morphological properties, 
globally-integrated spectra, and line ratio diagnostics of 
these galaxies as derived from the integral-field data. In 
§5 we use extant broadband photometry and UV spec- 
troscopy of the target galaxies to model the stellar popu- 
lations and baryonic masses of the galaxies. We present 
spatially resolved maps of the ionized gas kinematics in 
§6, discussing the velocity field of each of our targets in 
the context of theoretical models and the physical prop- 
erties given in §3-5. Finally, in §7 we compare our re- 
sults to those of Forster Schreiber et al. (2006) and 
other groups and discuss what general conclusions can 
be drawn about galaxy formation in the early universe. 

We assume a standard ACDM cosmology in which 
iJo = 71 km s-i Mpc-i, 0.27, and f^A^^ 0.73. 

2. OBSERVING AND DATA REDUCTION 

2.1. Sample Selection and Observational Strategy 

Observations were performed using the OSIRIS (OH 
Suppressing InfraRcd Imaging Spectrograph; Larkin et 
al. 2006) integral-field spectrograph in conjunction with 



the Keck II LGSAO system (Wizinowich et al. 2006, van 
Dam et al. 2006). As described in detail by Larkin et 
al. (2006), OSIRIS is a lenslet-type spectrograph with a 
2048 x 2048 Hawaii II detector array, spectral resolution 
R 3600, and a complement of reimaging optics which 
permit sampling of the field with spectral pixels (spaxels) 
ranging in angular size from 20 — 100 milliarcseconds 
(mas). 

We select star forming galaxies at redshifts z ~ 2 — 3 
(for which Ha and [O III] A5007, hereafter simply [O ill] , 
respectively fall into the observed-frame K-hand) from 
the spectroscopic redshift catalog of Steidel et al. (2003, 
2004). Based on the IFU simulations of Law et al. (2006) 
which explore the flux requirements necessary for rea- 
sonable detection within a few hours of integration using 
OSIRIS, we focus on those galaxies which are known to 
have either nebular line fluxes greater than ^ 10^^^ erg 
cm~^ based on previous long-slit observations (Erb et 
al. 2003, 2006a) or high star formation rates (SFR) as es- 
timated from rest-frame UV photometry (for one galaxy 
[Q0449-BX93] not targeted with long-slit spectroscopy). 
Physical constraints also dictate that we preferentially se- 
lect galaxies which are within 60 arcseconds of a suitably 
bright (i? ^ 17 mag) tip-tilt reference star (as required 
by the LGSAO system), and which lie at spectroscopic 
redshifts such that nebular emission lines fall between the 
strong spectroscopically unresolved night-sky OH emis- 
sion features which dominate the near-IR background. 

At if-band wavelengths between the OH-emission fea- 
tures, thermal radiation from warm optical components 
in the light path is the dominant contributor to the total 
background flux (see Law et al. 2006). Since successful 
observation of our target galaxies requires performance 
at the limit of OSIRIS' capability, we use the second- 
largest spaxel scale^ of 50 mas lenslet"-'^, for which the 
balance between thermal background and detector noise 
is optimal for detection of extremely faint objects. Us- 
ing narrow-band order sorting filters (AA 100 nm), 
the OSIRIS field of view in this configuration is approx- 
imately 2x3 arcseconds. 

Observations were performed on UT 2006 June 5 and 
2006 October 4-5; the June data were obtained under ex- 
ceptional photometric conditions with visual-band seeing 
estimates varying between 0.4 and 0.6" while the October 
data were obtained in moderate conditions with seeing 
varying between 0.7 and 1.0" over the duration of the 
observations. In Table [T] wc list the three targets dis- 
cussed herein as part of our pilot program. Each target 
was observed with the aid of the Keck LGSAO system, 
typically using a magnitude i? ~ 15 — 16 tip-tilt (TT) 
reference star with an angular separation from the sci- 
ence target of less than an arcminute for a correction 
strehl of greater than 30% at A = 2/im and an approxi- 
mately circular PSF with FWHM ~ 75/125 mas for the 
June/October data respectively. ^ 

A typical observing sequence was established as fol- 

^ Note that 50 mas is approximately the Keck difTraction limit 
at A = 2/im. 

* The 75 mas measurement is strictly an upper limit since this 
PSF is undersampled using 50 mas spaxels. Since science obser- 
vations are performed with the TT star off-axis the PSF of the 
science observations will not perfectly match that of the TT star; 
however, this generally does not greatly increase the effective PSF 
(see Liu et al. 2006). 
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lows. First, we obtained short (~ 1 minute) observations 
of the TT reference star for a given science target to refine 
the pointing model, provide an estimate of the PSF, and 
serve as an approximate check for later flux calibration. 
The telescope was then blind-offset to the science target 
using relative offsets calculated from deep ground based 
optical imaging of the target fields. The position angle of 
OSIRIS for these observations was chosen to ensure that 
the tip-tilt reference star fell within the unvignetted field 
of view of the LGSAO system. Each target was observed 
in sets of two 15 minute^ integrations with the target 
galaxy alternately placed 0.7 arcseconds above and below 
the center of the IFU (positions A and B). Small dithers 
were introduced around these positions on each succes- 
sive set of exposures. This 30 minute exposure sequence 
was repeated for ^ 2-4 hours as necessary to achieve a 
high-quality detection; typically it was possible to verify 
detection of target galaxies in the difference of two 15 
minute exposures. 

2.2. Reducing IFU Data 

Data reduction was performed using a multi-step pro- 
cess consisting of both Keck/OSIRIS data reduction 
pipeline (DRP) routines and custom IDL scripts tailored 
for use with faint emission line spectra. Since OSIRIS is 
a relatively new scientific instrument and standard data 
reduction methods have yet to be generally accepted, we 
present our current data reduction method in detail. 

We divide the exposures into two sets taken at posi- 
tions A and B, and use DRP routines to median combine 
each set to produce a sky frame. These sky frames are 
then subtracted from each of the respective object frames 
(i.e. the sky produced by stacking observations at posi- 
tion B is subtracted from the observations at position A 
and vice versa). After sky subtraction, individual expo- 
sures are extracted from the raw two-dimensional format 
into three-dimensional data cubes. This DRP extraction 
algorithm is complicated and non-intuitive since spectra 
from individual spaxels physically overlap on the detec- 
tor; we refer the reader to the OSIRIS user manual for 
further information. Due to the highly time-variable na- 
ture of the IR background the residual flux pedestal in 
these sky-subtracted data cubes is often unsatisfactorily 
high and we therefore perform a second-pass background 
subtraction on each of these cubes as follows. Using cus- 
tom IDL routines we calculate the median pixel value 
in each spectral channel and subtract this value from all 
pixels within the channel to ensure a zero-flux median 
in all spectral slices throughout the data cube. This ef- 
fectively removes flux pedestals resulting from variable 
sky brightness, although small variations in the back- 
ground persist across the spatial dimensions of each cube. 
As OSIRIS data reduction algorithms are progressively 
refined it may prove possible to further suppress back- 
ground fluctuations, achieving effectively deeper integra- 
tions. 

These sky-subtracted data cubes are spatially regis- 
tered using the telescope steering system (TSS) offset 
coordinates contained within the image headers and com- 
bined using a 3ct clipped mean algorithm. The resulting 

® Due to the read noise characteristics of OSIRIS, long individual 
exposures are desirable in order to optimize the S/N ratio of the 
observations. 



cube is spatially oversampled by a factor of two in each 
dimension (i.e. mapped to pixels 25 mas in size), rotated 
with flux-conserving resampling to a standard (North up. 
East left) orientation, and the image slice at each spec- 
tral channel convolved with a spatial gaussian kernel of 
FWHM 3.2 oversampled pixels (80 mas, giving an ef- 
fective PSF FWHM of 110/150 mas for June/October 
data respectively) in order to increase the S/N ratio of 
each spectrum. Wavelength calibration is automatically 
performed by the OSIRIS DRP; we determine that this 
calibration is generally accurate to within ~ 0.1 A based 
on fits to the centroids of strong OH sky lines with re- 
gard to fiducial vacuum wavelengths given by Rousselot 
et al. (2000). Finally, all spectra are converted to the 
heliocentric rest frame using the relative motion solu- 
tion given by the IRAF / rv correct package. Based on the 
measured FWHM of the (unresolved) OH lines, we con- 
clude that the effective spectral resolution of our data is 
R - 3600.1" 

Since we do not dedicate time in each observing se- 
quence to obtaining a clean "sky" observation, nega- 
tive residuals of the target galaxy are located about 1.4" 
above and below the positive image in the final stacked 
data cube, effectively halving our useful field of view. 
However, given the small size of the target galaxies (UV 
morphologies indicate that 90% are less than 1" in ex- 
tent; Law et al. 2007) we find this approach preferable 
for maximization of on-source integration time by elim- 
inating the need to frequently nod off-target in order to 
adequately characterize the background flux (although 
c.f. discussion by Davies 2007). 

3. NEBULAR MORPHOLOGIES 

We produce maps of the ionized gas morphology of 
the three target galaxies by collapsing each data cube 
about the 4-5 spectral channels corresponding to the 
peak emission wavelength of the spectral line (Ha for 
Q1623-BX453 and Q0449-BX93, [Oiii] for DSF2237a- 
C2). As illustrated in Figure [T] (top left panel), Q1623- 
BX453 consists of an Ha bright nucleation^^ with a lower 
surface brightness region extending towards the south- 
west and an extremely faint feature on the western edge. 
This faint feature is only barely detected and may be 
an artifact of imperfect sky subtraction. While fairly 
centrally concentrated, the bright emission region is in- 
consistent with a point source, having a diameter of 200 
mas (1.7 kpc) at the half peak- flux contour compared 
to (an upper limit of) 110 mas for the PSF of the tip- 
tih reference star. Q0449-BX93 (Fig. [H middle left 
panel) has a morphology quite similar to that of Q1623- 
BX453, consisting of a bright nucleated region with low 
surface brightness emission extending towards the north- 
west. Again, the bright region is inconsistent with a 
point source, having diameter 180 mas (1.5 kpc) at the 
half peak-flux contour. In contrast, DSF2237a-C2 (Fig. 
[11 bottom left panel) has a more elongated morphology 
than either Q1623-BX453 or Q0449-BX93, with flux dis- 

This represents the narrowest OH line observed; many ap- 
pear to be broader because the OH forest is typically composed of 
unresolved doublets separated by less than ~ lA. 

We adopt the term "nucleation" to qualitatively describe a 
concentrated region of high surface brightness that might naively 
be described as the "nucleus" of a given galaxy (see discussion in 
Law et al. 2007). 
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TABLE 1 
General Information 



Name RA {J2000) DEC {J2000) Observed Time=' Rtt^^ 9tt " ^seeing fpSF 

Q1623-BX453 16:25:50.854 +26:49:31.28 Jun '06 2.5 15.1 32 0.4 75/110 

Q0449-BX93 04:52:15.417 -16:40:56.88 Oct '06 4.5 15.8 51 0.6 125/150 

DSF2237a-C2 22:40:08.298 +11:49:04.89 Jun '06 1.5 16.0 24 0.4 75/110 



^ Total observing time in hours. R-hand magnitude of TT star for LGSAO correction.'^ Angular separation 
(arcseconds) of TT star from target.'^ Average ii'-band seeing (arcseconds) during observation.^ FWHM of the 
i<"-band PSF (mas) during on-axis TT star observation (before/after spatial smoothing respectively). 



tributed fairly evenly along the major axis. 

Given the elongated morphology of DSF2237a-C2, one 
traditional assumption would be that emission in this 
galaxy is predominantly from an inclined disk of ionized 
gas, and that the apparent ellipticity is due to perspec- 
tive foreshortening of a circularly symmetric disk. Un- 
der such an assumption, we fit a 2D gaussian to the Ha 
surface brightness map and conclude that with a ma- 
jor/minor axis ratio of rmaj/^mln — 

1.53, such a disk has 
a position angle oi 9 = 38° East of North and an incli- 
nation oi i = 49° away from the normal. 

While it is tempting to draw such physical inferences 
from the ionized-gas morphologies, these inferences are 
not necessarily robust given the irregularity of typical 
star-forming galaxies at redshift z ~ 2 — 3, which typ- 
ically consist of 1-3 nucleations connected by regions of 
low surface-brightness emission (as viewed in the rest- 
frame UV; see Fig. 2 of Law et al. 2007). Since the UV 
and ionized gas morphologies are strongly correlated (at 
least in the local universe; Gordon et al. 2004, Kenni- 
cutt et al. 2004, Lee et al. 2004, although c.f. Conselice 
et al. 2000) and trace the brightest regions of star for- 
mation unobscured by dust they are therefore a complex 
product of the strength and distribution of star form- 
ing regions, dust, viewing angle, and possible merger- 
induced irregularities, rendering them unreliable indica- 
tors of the galactic potential. Indeed, recent studies have 
found little to no correlation between morphology and 
major-axis velocity shear (Erb et al. 2004) or any of 
an assortment of other parameters such as stellar mass, 
star formation rate, or the strength of galactic-scale out- 
flows (Law et al. 2007). Even if the rotating gas disk 
hypothesis were correct, it is not obvious that estimates 
of the disk size and/or orientation based on the length 
of the major/minor axes are meaningful since the ob- 
served morphology of such a disk will be dominated by 
the clumpy distribution of star formation (see, e.g., theo- 
retical models by Noguchi 1999 and Immeli et al. 2004b). 

We quantify an average linear dimension of the nebu- 
lar line emission region by adopting a segmentation map 
which includes all spaxels (i.e. lenslets) for which nebular 
emission can be fit with a gaussian profile with S/N ratio 
(SNR) ~ 6, and which have a velocity of less than 300 km 
relative to systemic (see §6 for further explanation of 
these selection criteria). The effective size of this region 
is determined as follows. Assuming that the N elements 
in the segmentation map are distributed in an approxi- 
mately circular arrangement, the radius of this circle is 
given simply by r' — y^N/ir. The corrected radius r is 
then determined by subtracting off the radius of the stel- 
lar PSF in quadrature (Table [U this is typically a small 
correction of less than ^10%) and converting to physi- 



TABLE 2 
OSIRIS Morphologies. 



Galaxy 


















Q1623-BX453 


9.5 


± 


0.3 


1.74 


± 


0.03 


0.33 


1.48 


Q0449-BX93 


4.7 


± 


0.6 


1.23 


± 


0.07 


0.25 


1.32 


DSF2237a-C2 


3.3 


± 


0.3 


1.03 


± 


0.04 


0.29 


3.33 



^ Area of nebular emission (kpc^). Uncertainty rep- 
resents half the PSF correction.'^ Radius of nebular 
emission (kpc). Uncertainty represents half the PSF 
correction.*^ Gini.*^ Multiplicity. 



cal kpc at the redshift of the source. The PSF-corrected 
emission area is defined as / = 7rr^. Adopting the non- 
parametric gini- multiplicity (G — 5*) classification scheme 
described by Law et al. (2007), we note that the ionized- 
gas morphologies of our three target galaxies (Tabic ^ 
closely resemble the typical rest-frame UV morphologies 
of z ^ 2 — 3 star forming galaxies in the GOODS-N (Law 
et al. 2007) and Q1700 (Peters et al. 2007) fields as seen 
by HST-ACS, albeit with limiting star formation rate 
density roughly five times higher than that of the HST 
data. 

4. GLOBAL SPECTRA 

4.1. Flux Calibration 

Integrated spectra (shown in Fig. ^ were generated for 
each source by convolving the data cube during reduction 
with a broad 150 mas spatial gaussian kernel, which con- 
siderably improves the S/N ratio (SNR) of each spaxel 
at the cost of suppressing spatial information. All spax- 
els with SNR > 6 were combined to produce a "clean" 
integrated spectrum, while all spaxels within a radius of 
about 2r were combined to produce a "total fiux" inte- 
grated spectrum. The first of these is of higher quality 
and is used to determine the relative flux in individual 
emission lines, while the second spectrum better samples 
the faint flux distribution and is used to flx the abso- 
lute fiux calibration. Wc eliminate any residual gradients 
as a function of wavelength by subtracting off a first- 
order spline flt to the integrated spectrum using stan- 
dard IRAF techniques (with iterative rejection to exclude 
emission lines from the spline fit). We estimate the flux 
calibration of our integrated emission-line spectrum us- 
ing detailed modelling of the OSIRIS + Keck LGSAO 
system -I- atmospheric throughput as a function of wave- 
length based on the simulations of Law et al. (2006) and 
refined to accomodate recent measurements of system 
throughput parameters. The results of this wavelength- 
dependent calibration method agree to within ~ 25% 
with those obtained using the broadband photometric 
zero-point calibrations published in the OSIRIS manual 
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Fig. 1.— OSIRIS maps of (left to right) nebular emission line flux density (Ha for Q1623-BX453 and Q0449-BX93, [O III] for DSF2237a- 
C2), velocity, and velocity dispersion for the three target galaxies. Stretch on plots of emission line flux density is linear and runs from zero 
to the peak flux density of each galaxy, linear contours in flux density are overplotted on the velocity and dispersion maps. The FWHM 
of the PSF after smoothing is 110 mas in Q1623-BX453 and DSF2237a-C2 and 150 mas in Q0449-BX93, indicated by solid lines in the 
left-hand panels. Individual pixels measure 25 mas, the total field of view (0.75 X 0.75 arcseconds) is comparable to the seeing-limited 
spatial resolution of previous work. This field of view corresponds to 6.3, 6.4, and 5.7 physical kpc at the redshift of Q1623-BX453, 
Q0449-BX93, and DSF2237a-C2 respectively 

quantities for these various parameters roughly represen- 
tative of our observing conditions, we expect to recover 
~ 60% of the light from a given target source. This per- 
centage is roughly consistent with that estimated using 
measurements of the average flux density in narrow-band 
observations of our TT stars compared to Kg magnitudes 
drawn from the 2MASS database. 

Given all of the uncertainties inherent in this boot- 
strapped flux calibration, we estimate that the system- 
atic flux uncertainty for a given source is ~ 30%. In con- 
trast, the relative flux uncertainty between two emission 
lines (i.e. the uncertainty in the line flux ratio) is con- 
siderably lower and dominated simply by the root-mean 
square (RMS) deviation of the integrated spectrum. For- 
mally, we adopt as the la uncertainty in each emission 
line flux the error due to RMS fluctuations per spec- 



(roughly as close an agreement as may reasonably be ex- 
pected since these zero-points average the transmission 
function throughout the bandpass). The greatest uncer- 
tainty in flux calibration arises however in the additional 
requirement that we correct for the unknown amount of 
flux lost to the seeing halo of the AO-corrected PSF. 
For a point source with a 35% strehl correction approx- 
imately 35% of the total light is contained within the 
diffraction-limited core of FWHM ~ 50 mas while the re- 
maining 65% is spread throughout a seeing-limited halo. 
The percentage of total flux measured for a given obser- 
vation is therefore a complex function of the strehl of the 
AO correction (which depends in turn upon the magni- 
tude/separation of the TT star and the wavelength of 
observation) , the quality of the seeing, and the physical 
structure of the source. Based on simulations adopting 
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tral channel compounded over the measured width of the 
emission hue. Similar estimates of the RMS uncertainty 
may be used to place Scr upper limits on the strength of 
undetected emission lines using a reasonable guess at the 
likely intrinsic width of the undetected line. 

Based on these calibrations and observations of the 
three target galaxies we determine that /ii,„ = 3 x 10~^^ 
erg s~^ cm~^ arcsec"^ represents the 6(7 Ha surface 
brightness threshhold for our observations (i.e. the 
faintest believable feature for two hours of integration 
at A ~ 2/x). At redshift z — 2.5 this is equivalent to a de- 
tection threshhold of SFRum = IM© yr~^ kpc~^ (adopt- 
ing the Ha flux to SFR calibration described in §5.2).^^ 
Given that a star formation rate density of ^ 0.1 Mq 
yr~^ kpc~^ represents the threshhold between starburst 
galaxies and the centers of normal galactic disks in the 
local universe (Kennicutt 1998), we do not expect to de- 
tect local-type disk galaxies in the course of our study 
even if such galaxies existed at high redshift. 

In Table [3] we compile the systemic, UV absorption- 
line, and Lya emission redshifts and nebular line fluxes 
of the target galaxies. In both cases for which previous 
NIRSPEC observations were obtained (Q1623-BX453 
and DSF2237a-C2) the systemic redshifts agree to within 
Az — 0.0008 or about 50 km s""'^ (i.e. consistent to 
within roughly one fifth the velocity resolution of NIR- 
SPEC [i? - 1400]). The integrated Ha flux of Q1623- 
BX453 is 1.7 times greater in the OSIRIS data than the 
NIRSPEC, consistent with the approximately 50% slit- 
loss correction factor estimated by Erb et al. (2006a). 
However, the [O ill] flux of DSF2237a-C2 appears to be 
a factor of ~ 2 weaker in the OSIRIS data than the NIR- 
SPEC. The first, and perhaps most likely, explanation 
for this discrepancy is that the NIRSPEC fiux observed 
is strongly sensitive to differential slit losses between the 
standard star and target galaxy observations, which can 
drive flux errors in both directions depending on atmo- 
spheric conditions and the quality of the background sub- 
traction in each case. However, it is also possible that 
this galaxy contains a large fraction of its total flux in 
the form of extended low surface brightness emission un- 
detected with OSIRIS but detected in the more sensitive 
long-slit observations. 

4.2. Global Chemistry 

In the case of Q1623-BX453, both [Nii] A6549 and 
A6584 (the stronger A6584 line is hereafter simply [Nil]) 
are detected reliably in the integrated galaxy spectrum 
(see Figure[2|), allowing us to estimate the average metal- 
licity of the galaxy. Defining the flux ratio N2 = 
log([Nli]/Ha), we follow Pettini & Pagel (2004) who es- 
timate the oxygen abundance as 

12 -I- log(0/i/) = 8.90 + 0.57 X 7V2 (1) 

Adopting the [N ii] flux tabulated in Table [4] we find 
N2 = -0.48 ± 0.03, and therefore 12 + log(0/H) = 
8.63 ± 0.18 (this uncertainty is dominated by the 0.18 
dex systematic uncertainty in the calibration given in 

Note that this is the threshhold for detection, the threshhold 
for deriving reliable kinematic information is approximately 4 Mq 
yr^^ kpc~-^. For comparison, the limiting star formation rate den- 
sity in the rest-UV morphological analysis presented by Law et al. 
(2007) is ~ 0.2 Mq yr'^ kpc^^. 



Equation [T]) , consistent with near-solar enrichment (12 
+ log(0/H)o = 8.66; Asplund et al. 2004). This is 
consistent with the results of Shapley et al. (2004), 
who find 12 + log(0/H) = 8.60 ± 0.18 based on long- 
slit spectroscopy. As noted previously by Shapley et al. 
(2004) however, the FWHM of the [Nil] emission line is 
roughly 30% larger than that of Ha, which can be a sig- 
nature of shock ionization (Lehnert & Heckman 1996b). 
Since Q1623-BX453 shows no evidence of AGN activity 
in either rest-frame UV spectra (§4.3) or long-wavelength 
deviations from simple stellar population models (§5.1), 
and the measured value of A''2 is consistent with that 
observed for normal star-forming galaxies (Shapley et al. 
2004; Erb et al. 2006c) we see no concrete evidence that 
this high ratio is due to shock heating or physical pro- 
cesses other than simple enrichment. Erb et al. (2006c) 
have considered this question in greater detail using a 
[Oiii]/H/3 vs. [Nil] /Ha diagnostic diagram and find that 
while the [O iii]/H/3 ratio is greater for a given [Nil] /Ha 
ratio in Q1623-BX453 than in local galaxies, this offset 
may be explained by a harder ionizing spectrum or a 
higher electron density (see also discussion by Shapley et 
al. 2005b) and does not provide concrete evidence for a 
central AGN. 

Although the spectra are too faint to determine the 
precise [N ii] morphology, we note that [N ii] appears to 
be coincident with the Ha peak and varies slightly in rela- 
tive strength across the face of the galaxy. Figure [3] plots 
individual spectra of the northeastern and southwestern 
quadrants and illustrates the stronger [N ii] lines present 
in the northeastern region. Taking careful measurements 
of the respective spectra we find (Table [5]) that the fiux 
ratio of the bright core is N2 = —0.46 ± 0.03 (consistent 
with the overall spectrum; as expected since it domi- 
nates the total fiux of the source) while the flux ratio of 
the fainter region is A^2 = -0.67 ± 0.09. This ^ 2.3cr 
difference (~ 0.1 dex in log(0/H))offers mild evidence 
to suggest that the distribution of heavy elements in 
Q1623-BX453 is not perfectly well-mixed, suggesting ei- 
ther that heavy elements from the latest burst of star for- 
mation (produced primarily in the northeastern region) 
have not yet had time to be fully distributed throughout 
the galaxy, or that the two regions may have had dis- 
tinct star formation histories (as might be expected for 
instance in a merger-type scenario). 

Since Q0449-BX93 is considerably fainter than Q1623- 
BX453 it is not possible to detect [Nil] emission confi- 
dently. We place a 3a detection limit on /jj^ jjj (Table[3]) 

and conclude that an upper bound on the metallicity is 12 
+ log(0/H) < 8.37 (i.e. less than 1/2 solar). Given the 
monotonic relation between stellar mass and metallicity 
found by Erb et al. (2006c) for similar galaxies at z ~ 2, 
this metallicity limit implies that M^, < 10^*^ Mq. Unfor- 
tunately, it is not possible to constrain the metallicity of 
DSF2237a-C2 as we have only observed the wavelength 
interval surrounding [O ill] A5007A. 

4.3. Outflow Properties: UV Spectroscopy 

In Figure m we plot the Keck-LRIS UV spectra of the 
three target galaxies shifted into the rest-frame as de- 
fined by the OSIRIS nebular redshifts. The spectrum 
of DSF2237a-C2 was obtained as part of a survey of 
z ~ 3 LBGs (Steidel et al 2003; see also Giallongo et 
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Fig. 2.— OSIRIS spectra of Q1623-BX453, Q0449-BX93, and DSF2237a-C2 collapsed over the spatial extent of each galaxy. Vertical 
dotted lines denote the wavelengths of [N ii] A6550, Ua, and [N ii] A6585 emission (for Q1623-BX453 and Q0449-BX 93) and [O in] A4960 
and A5008 emission (for DSF2237a-C2). Lines for [N ii] and [O iii] A4960 are drawn at the fiducial wavelengths expected based on the Ha 
and [O III] A5008 emission redshifts. 



TABLE 3 

Primary Spectroscopic; Characteristics 



Galaxy 




Z 

^-■ncb 


^abs 








^[OIII]'' 










Q1623-BX453 


20888.30 


2.1820 


2.1724 


2.1838 


23.1 ± 0.5 


18.0 ±0.4 






92 ± 8 




17 ± 2 


Q0449-BX93 


19737.80 


2.0067 


2.004 




9.3 ±0.3 


4.1 ±0.1 






72 ± 12 




7 ± 2 


DSF2237a-C2 


21621.80 


3.3172 


3.319 


3.333 






10.1 ± 0.5 


2.2 ±0.1 


101/79 ± 15 


70 ± 10 


5 ± 1 



^ Vacuum wavelength of peak emission in Angstroms.*^ Vacuum redshift of nebular emission line: Ha for Q1623-BX453 and Q0449-BX93, [O in] A5007 
for DSF2237a-C2.*^ Vacuum heliocentric redshift of rest-frame UV interstellar absorption lines. Vacuum heliocentric redshift of Lya emission.*^ Emission 
line flux in units of 10""'"'^ erg s""*" cm~^. Uncertainties quoted arc based on random errors, global systematic uncertainty is ~ 30%.^ Extinction-corrected 
emission line luminosity in units of 10^^ erg s""'". Uncertainties quoted are based on random errors, global systematic uncertainty is ~ 30%.^ Velocity 
dispersion in units of km s~^. Two values are given for DSF2237a-C2, the overall dispersion of the integrated spectrum and the median velocity dispersion 
in a given spaxeh^ Circular velocity in units of km s~^/ Dynamical mass in units of IO^Mq. 
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Fk;. 3. OSIRIS spectra of the northeastern (NE) and southwestern (SW) regions of Q1623-BX453. Vertical dotted lines denote the 

wavelengths of [N ii] A6550, Ha, and [N ii] A6585 emission. 
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TABLE 4 

Secondary Spectroscopic Characteristics. 



Galaxy F[o m] (A4960)^ F[Nn] (A6549)^ F[n„](A6585)^ 

Q1623-BX453 ... 3.3 ±0.4 7.7 ±0.5 

Q0449-BX93 ... < 1.1 < 1.1 

DSF2237a-C2 2.5 ± 0.5 



^ Emission line flux in units of 10 erg s ^ cm ^. Limits arc 3cr limits. 
Uncertainties quoted are stochastic, global systematic uncertainty is ~ 
30%. 



TABLE 5 

Q1623-BX453 Regional Emission Line Characteristics''. 



Region 




/[NII]A6549 


/[NII]A6585 


N2^ 


NE 


14.8 ±0.4 


1.9 ± 0.4 


5.1 ±0.4 


-0.46 ±0.03 


SW 


6.6 ±0.4 


< 1.0 


1.4 ±0.4 


-0.67 ±0.09 



^ Emission line flux in units of 10 erg s ^ cm ^. Limits arc So" 
limits.'' N2 = log([Nn]/HQ). 



al 2002), while Q1623-BX453 and Q0449-BX93 were ob- 
tained more recently using LRIS-B, described by Steidel 
et al (2004). Q1623-BX453 was included on a slit mask 
receiving a total of 5.5 hours integration time, whereas 
the other spectra were obtained using shorter, survey- 
mode integration times of 1.5-2 hours. 

The spectra of all three galaxies are broadly simi- 
lar to high-quality composite spectra of the 2^2 — 3 
star- forming galaxy sample (e.g. Shapley et al. 2003); 
these spectra exhibit strong interstellar absorption fea- 
tures (e.g. Si II, O I, C II, Si iv, C iv) and Lyman a 
that may appear in either emission or absorption. It 
is typical of star forming galaxies at these redshifts for 
the interstellar lines to be blue-shifted, and Lya emis- 
sion red-shifted (when present), by several hundred km 
with respect to the systemic redshift; these features 
are generally interpreted as scattered emission from the 
far-side and absorption in the near-side of an expanding 
shell of outflowing gas (see, e.g., Steidel et al. 2007). The 
mean values for a sample of ~ 100 galaxies in the red- 
shift range 1.9 < z < 2.6 are Awabs = —164 km s~"'^and 
Afcm = +445 km s~^ for the centroid of the interstellar 
lines and Lya emission relative to Ha (which tends to 
trace gas at the systemic redshift), respectively (Steidel 
et al. 2007). 

The spectra of Q1623-BX453 and of Q0449-BX93 are 
typical with regard to the sense of the velocity offsets 
observed, as summarized in Table 3. The centroids of 
the interstellar lines of Q0449-BX93 are blue-shifted by 
~ —270 km s~^ with respect to the Ha redshift, placing 
it within the top 10 — 15% in terms of this quantity, 
and making it very similar to the well-studied spectrum 
of cB58 (Pettini et al 2000, 2002). We note in pass- 
ing that relatively strong Lya, CiV, and Si IV absorption 
lines at Zabs = 2.0074 are detected in the Keck/HIRES 
spectrum of the Zem = 2.677 QSO Q0449-1645. The 
difference from the systemic redshift of Q0449-BX93 is 
only 70 km s""'^, even though the projected separation 
between the QSO and galaxy sightlines is 44" , or 373 
(proper) kpc. 

The far-UV spectrum of Q1623-BX453 is more remark- 
able, in that the blue-shifted interstellar lines have a cen- 



troid offset of Awabs ^ -900 km s^^ , and a FWHM of 
~ 1200 km s~^ , both values being by far the largest ob- 
served in the sample of ~ 100 mentioned above. The 
Lya emission redshift exhibits a relatively modest shift 
of Avcm = +170 km s-i . The spectrum of Q1623-BX453 
also has unusually strong high-ionization interstellar ab- 
sorption features, including NV AA1238,1242, and many 
of the absorption features have complex velocity struc- 
ture, with two dominant components at —600 and —1100 
km . These unusual spectral features might lead one 
to suspect an AGN is involved, but there are no sig- 
nificant emission lines aside from Lya in the high qual- 
ity spectrum, and (as discussed further below) there is 
no evidence from the full spectral energy distribution to 
suggest the presence of anything but stellar energy. In 
spite of the rather violent kinematics of the outflowing 
interstellar gas, the Ha emission line kinematics show no 
hint of being affected, presumably because any Ha emis- 
sion associated with the outflowing material is of much 
lower surface brightness than what we have detected with 
OSIRIS. 

The spectrum of DSF2237a-C2 has typically strong in- 
terstellar absorption lines, but in this case close exam- 
ination reveals that they are actually redshifted slightly 
with respect to the nebular redshift, with Awabs — +100 
km s""'^ , while the Lya emission line is redshifted by an 
unusually large ~ 1100 km s~^ . It turns out that such 
kinematics, where the interstellar line centroids are con- 
sistent with zero velocity with respect to the nebular sys- 
temic velocity, are observed in ~ 15% of z ~ 2—3 galaxies 
(Steidel et al 2007), and it is quite typical for such objects 
to have unusually large Lya emission offsets. Of possible 
significance is that interstellar absorption line centroid 
velocities consistent with zero occur relatively frequently 
(5 out of 7) for galaxies in the Erb et al (2006c) sample 
having measured velocity shear from NIRSPEC Ha spec- 
tra. As discussed further by Steidel et al (2007), these 
observations suggest a connection between the presence 
of velocity shear (as in the case of DSF2237a-C2; see §6.3) 
and increased absorption close to zero systemic velocity, 
though at present it is not possible to distinguish be- 
tween orientation effects (e.g., outflows collimatcd along 
the rotation axis) and real differences in the dynamical 
structure of the galaxies. In any case, OSIRIS is much 
more sensitive to the measurement of velocity shear in 
the galaxies, and any correlation between the kinematics 
of gas associated with outflows and the kinematics of the 
H II regions should become clearer as the sample size 
increases. 

5. STELLAR POPULATIONS, STAR FORMATION RATES 
AND GAS MASSES 

5.1. Stellar Masses and Stellar Populations from SED 

Modeling 

As part of an ongoing survey program we have ob- 
tained broadband photometry (Table [S]) in each of our 
target fields. The wavelength coverage of these data 
varies from field to field: aU fields (Q1623, Q0449 and 
DSF2237a) have ground-based optical J7„G72. photome- 
try which serves as the basis for the optically-selected 
galaxy catalog (Steidel et al. 2003, 2004); Q1623 and 
DSF2237a have deep near-IR J and Kg photometry, de- 
scribed by Erb et al. (2006b) and Shapley et al. (2001) 
respectively; and Q1623 has also been imaged in the mid- 
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Fig. 4. — Flux calibrated UV spectra of the three target galaxies. Spectra have been shifted to the systemic rest frame using the nebular 
line redshifts given in Table |3] The identifications and wavelengths of the strongest spectral features are indicated by vertical dotted lines. 
Spectra of Q0449-BX93 and DSF2237a-C2 have been smoothed with a 3-pixel boxcar filter. 



IR by the Spitzer Space Telescope, with both IRAC at 3-8 
and MIPS at 24 ^m. 

Wavelength coverage from the rest-frame UV to near- 
IR (in Q1623) or UV to optical (DSF2237a) enables us to 
model the spectral energy distributions of Q1623-BX453 
and DSF2237a-C2 and obtain estimates of these galaxies' 
stellar masses, ages, reddening and star formation rates 
(the Q0449 data are restricted to the rest-frame UV and 
are insufficient for such modeling). Previous modeling 
of the SEDs of these two galaxies has been presented by 
Erb et al. (2GG6b) and Shapley et al. (2001) respectively; 
we redo the modeling here to add the Spitzer IRAC data 
in the case of Q1623-BX453 and update the method for 
consistency in the case of DSF2237a-C2. In both cases 
the Ks magnitudes have been corrected for line emission 
according to the values determined in §4.1 (the correc- 
tions are ~ 0.2 mag in both cases). 

The modeling procedure has been described in detail 
by Shapley et al. (2005a). Briefly, we use Bruzual & 
Chariot (2003) models with a constant star formation 



rate, solar metallicity, and a Chabrier (2003) IMF. The 
resulting stellar masses for Q1623-BX453 and DSF2237a- 
C2 are = 3.1 x lO^" Mq and A'h = 2.9 x lO^" Mq 
respectively. The remaining stellar population parame- 
ters are given in Table [6l and the best-fit spectral energy 
distributions are shown in Figure El 

5.2. Star Formation Rates 

A direct estimate of the instantaneous star formation 
rate may be derived from the observed Ha luminosity 
since it traces the ionizing radiation produced by young, 
massive stars. We adopt the calibration derived by Ken- 
nicutt et al. (1994; summarized in Kennicutt 1998) 
whereby 

L(Ha) 



SFR (Mq yr-^ - 



X 0.56. 



(2) 



1.26 X 1041 erg s"! 
where we have introduced the factor of 0.56 to convert 
to the Chabrier (2003) IMF. 

This enables direct calculation of the star formation 
rate from our observed Ha fluxes for Q1623-BX453 and 
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TABLE 6 

Photometric Data and SED Model Parameters 



Galaxy 






7^=' 








"14. Sm"" 










Age* 


E(B - VY 


Q1623-BX453 


24.85 


23.86 


23.38 


21.41 


19.76 [20.00] 


21.38 


21.37 


21.33 


21.92 


195 


3.1 


404 


0.245 


Q0449-BX93 


23.62 


23.17 


22.99 
















<lh 




0.125 


DSF2237a-C2 




24.68 


23.55 


22.08 


20.53 [20.71] 












2.9 


255 


0.175 



^ Un,G,7Z magnitudes are AB. J, Kg magnitudes are Vega. Values in brackets for Ks magnitudes represent continuum magnitudes corrected for 
line emission. Spitzer-IRAC magnitudes (AB standard).'^ Spitzer-MIPS 24fi flux (AiJy)-° Stellar mass (IO^^A/q) from the best-fitting constant star 
formation model, using the Chabricr (2003) IMF. Typical uncertainty < ctm^/M^ >— 0.4.^ Stellar population age (Myr) from SED fitting, typical 
uncertainty < fTAge/Agc >^ 0.5.S E{B - V) from the best-fit SED for Q1623-BX453 and DSF2237a-C2; E{B - V) is estimated for Q0449-BX93 from 
the G — "R, color, assuming an age of 570 Myr.^ Estimated from the metallicity limit (§4.2) using the mass-met allicity relation of Erb et al. (2006c). 



TABLE 7 

Star Formation Rates and Gas Masses 



Galaxy 




SFRhc" 


SFRuv*" 


SFRuv'' 




SFRsed'* 




SSFR^ 


Mgas*" 






Uncorrected 


Corrected 


Uncorrected 


Corrected 






Mean 


Peak 






Q1623-BX453 


36 


77 


8 


81 


96 


77 


8 


20 


2.4 


0.44 


Q0449-BX93 


12 


18 


12 


39 






4 


7 


0.7 


> 0.41 


DSF2237a-C2 


23 


49 


18 


95 




112 


15 


25 


1.3 


"b.32 



^ Star formation rate (-A/0 yr ^) estimated from Ha flux. Typical uncertainty < ctsfr/SFR >— 0.3.^ Star formation rate {Mq yr ^) 
estimated from rest-frame UV continuum flux.^ Star formation rate (Mq yr~^) estimated from rest-frame 8fi flux.'^ Star formation rate 
{Mq yr~^) derived from SED modeling. Typical uncertainty < ctsfr/SFR >— 0.6.^ Star formation rate surface density {Mq yr~^ 
kpc~^).^ Molecular gas mass (IO^^A/q). Typical uncertainty < (7m /M >— 0.2.^ Gas fraction fi — A/gas/(-A/gaB + A/*). Typical uncertainty 
< o-^//i >= 0.35. 




20 



U G R 



K 



22 



24 - 



Q1623-BX453 
z = 2.1820 
Mass = 3,lxlOio M„ 
Age=404 Myr 



2 4 6 

Observed Wavelength (//m) 



26 




DSF2237a-C2 
z = 3.3172 
Mass = 2. 9x1010 M„ 
Age = 255 Myr 



1 2 
Observed Wavelength (/xm) 



Fig. 5. — The best-fit constant star formation (CSF) model (solid black line) is overplotted against the observed spectral energy 
distribution for Q1623-BX453 and DSF2237a-C2. Colored points represent ground-based optical and near-IR photometry, black points are 
based on Spitzer-TRAC observations. The _fC-band photometry for both galaxies is shown corrected for nebular line emission (see discussion 
in §5.1). Values given for stellar mass and population age represent the values derived from the best-fit CSF model; typical uncertainties 
are given in Table [6] 

Q0449-BX93; we find SFRho = 36 Mq yr^^ for Q1623- 
BX453 and SFRho = 12 Mq yr"! for Q0449-BX93. Al- 
though Ha is redshifted beyond the X-band for z > 2.6, 
it is nonetheless possible to estimate a star formation 
rate for DSF2237a-C2 based on the observed [O III] flux. 
Using long-slit spectroscopy, Pettini et al. (2001) de- 
termined the H/3 flux to be 3.5 ± 0.4 x 10~^^ erg s~^ 
cm~^, and the ratio of the [O III] and H/3 fluxes to be 
Fm/Fioui] = 0.19 (their [O III]A5007 flux of 18 x 10 



17 



erg s^ cm^ is in rough agreement with the OSIRIS 
value of 10.1 ± 0.5 x 10~^^ erg s~^ cm~^, given the sig- 
nificant uncertainties in the calibration of long-slit spec- 
tra). While up falls outside the narrow wavelength 
range which we have observed with OSIRIS, we can use 
the ^H;3/-F[oiii] ratio observed by Pettini et al. (2001) to 
infer an H/3 flux of 2.0 x 10~^^ erg cm~^ from our 
observed [O III] value. We then assume a standard re- 
combination ratio of Fna/Fap = 2.75 (Osterbrock 1989) 
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to find FHa = 5.4 X IQ-^^ erg s^^ cm^^ ^nd SFRna = 23 
M0 yr~^. These SFRs are uncorrected for dust extinc- 
tion (which we discuss below) and have typical uncer- 
tainties ^ 30% (dominated by the uncertainty in the flux 
calibration described in §4.1). 

We use the values of E{B — V) determined by the 
SED modeling discussed in ijS.ll to correct the observed 
Ha flux of Q1623-BX453 and the inferred H/3 flux of 
DSF2237a-C2 for extinction using a modified Calzetti 
et al. (2000) starburst attenuation law as described 
by Erb et al. (2006b). For Q1623-BX453 we find 
E{B -V) = 0.245 for a corrected SFR of 77 M© yr^^, 
in agreement with the value of 77 Mq yr~^ we obtain 
from the SED modeling. For DSF2237a-C2 we use the 
best fit E{B -V) = 0.175 to correct the H/3 fiux in- 
ferred above, and then assume ^hq/-Fh/3 = 2.75 to find 
Fua = 11-3 X 10~" erg s'^ cm'^ and SFRsa = 49 
Mq yr^^. This is somewhat smaller than the value of 
112 M© yr~^ we find from the SED modeling, but not 
inconsistent with the factor of ~ 2 RMS found by Erb 
et al. (2006b) for SFRs determined from SED modehng 
and Ha emission, especially given the additional extrap- 
olations needed to calculate the nebular line SFR.^^ 

Because we lack the IR data to model the SED of 
Q0449-BX93, we estimate the extinction from the G — TZ 
color alone, assuming a constant star formation age of 
570 Myr, the median age^** of the z 2 sample as 
found by Erb et al. (2006a). Assuming this median 
age we find E{B — V) — 0.125 and a corrected SFR of 
18 Mq yr~^. The uncertainty in age has a negligible 
effect on the corrected SFR; for an age of 50 Myr, we 
find SFRhq — 20 Mq yr~^, while the maximum age of 
3.2 Gyr leads to SFRna = 17 M© yr"!. The corrected 
and uncorrected SFRs for the three galaxies are summa- 
rized in Table [71 and for reference we note that the mean 
extinction-corrected SFR of the z ^ 2 sample, from both 
Ha and the UV continuum, is ~ 30 M© yr^^ (Erb et al. 
2006a). 

By way of comparison, we also include in Table [7] esti- 
mates of the star formation rate obtained using UV and 
mid-IR photometry. In the former case, we use the rest- 
frame monochromatic luminosity at 1500A(L^; based on 
linear interpolation of the broadband photometry given 
in Table and the Kennicutt (1998b) conversion be- 
tween Li, and SFR: 



SFR (M© yr-i) 



7.14 X 1027 erg s-i Hz" 



X 0.56 



(3) 

where the factor of 0.56 has again been introduced to 
convert to the Chabrier (2003) IMF. After correcting for 
extinction using the Calzetti et al. (2000) starburst at- 
tenuation law, we find SFRuv = 81, 39, 95 M© yr~^ 
for Q1623-BX453, Q0449-BX93, and DSF2237a-C2 re- 
spectively. These values are in good agreement with 

Curiously, the original H/3 flux measured by Pettini et al. 
(2001) gives a better match to SED star formation rate than the 
extrapolated flux measured using OSIRIS, highlighting the possi- 
bility that OSIRIS may sometimes 'miss' low surface brightness 
flux. 

The values of E{B - V) allowed by the G - 7^ color are 
relatively insensitive to the age of the galaxy; for reasonable ages 
ranging from 50 Myr (approximately the dynamical timescale) to 
3.2 Gyr (the age of the universe at the redshift of the galaxy), 
E{B ~ V) varies from 0.164 to 0.113. 



SFRho for Q1623-BX453; the factor of - 2 difference 
between UV and Ha estimates for DSF2237a-C2 and 
Q0449-BX93 may be due (in the former case) to the un- 
certainties introduced in our conversion from [O ill] fiux 
to Ha and suggests (in the latter case) that Q0449-BX93 
may be a particularly young, blue, dust-free galaxy for 
which the conversion from (G — TZ) color to E{B — V) is 
imperfect. Using Spitzer-MlPS 24/i photometry, it is also 
possible to estimate the SFR for Q1623-BX453 using the 
rest-frame 8/i fiux given in Table [51 Adopting the cali- 
brations discussed in Reddy et al. (2006), we estimate 
SFRiR = 96 Mq yr-i under the Chabrier (2003) IMF, in 
good agreement with the UV, SED, and Ha estimates. 

5.3. Gas Masses and Gas Fractions 

We can also use the Ha luminosities in combination 
with the galaxy sizes measured in §3 to estimate the 
galaxies' gas masses and gas fractions. Such indirect 
techniques are required because emission from molecular 
gas is very difficult to detect in all but the most lumi- 
nous high redshift galaxies with current technology. We 
follow the procedure described by Erb et al. (2006a), 
benefitting here from the higher spatial resolution of the 
OSIRIS data, which provides better constraints on the 
galaxies' sizes. As described by Erb et al. (2006a), we 
use the empirical global Schmidt law (Kennicutt 1998), 
which relates the surface density of star formation to the 
surface density of molecular gas (in the starburst regime 
in which our galaxies fall; more quiescent star-forming 
galaxies at low redshift also have significant amounts of 
atomic gas). 

We determine the surface density of Ha luminosity 
(Shq) by dividing the extinction-corrected observed (or 
inferred, in the case of DSF2237a-C2) Ha luminosity 
by the projected linear area / of the galaxy as deter- 
mined in §3. The gas surface density Sgas is then given 
by the global Schmidt law, and the total gas mass by 
-^^gas — Sgas/. Gas masses inferred in this way are sub- 
ject to significant uncertainties beyond those associated 
with the measurements of fiuxes and sizes and the ~ 0.3 
dex uncertainty introduced by the scatter in the Schmidt 
law itself. It is not yet known whether the Schmidt law in 
this form holds at high redshift, although the one galaxy 
for which it has been tested, the gravationally lensed Ly- 
man break galaxy MS1512-cB58, appears to be consis- 
tent with the local form (Baker et al. 2004). It is also 
likely that the Schmidt law, which applies to only the 
cold gas associated with star formation, does not account 
for all of the gas; in a typical local disk galaxy, ~ 40% 
of the total gas mass is not included by the Schmidt law 
(Martin & Kennicutt 2001). 

The measured areas of Q1623-BX453, Q0449-BX93, 
and DSF2237a-C2 are / = 9.5 ± 0.3, 4.7 ± 0.6, and 3.3 ± 
0.3 kpc^ respectively. The size of Q1623-BX453 derived 
here (r — 1.74±0.03 kpc) is considerably smaller than the 
value of r = 4.4 ± 0.4 found by Erb et al. (2006b) from 
long-slit NIRSPEC observations, reflecting in part the 
slightly different definitions of r and the difficulty of esti- 
mating sizes from seeing-limited long-slit data. However, 
it is also likely that the luminous area of a given galaxy 
is genuinely greater in the Erb et al. (2006b) data since 
NIRSPEC is more sensitive than OSIRIS, and may thus 
detect lower surface-brightness features at larger radii. 
Galaxy size r and / as defined here are only meaningful 
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with regard to the detection threshhold of a particular 
study. 

Combining these size estimates with the extinction- 
corrected Ha luminosities, we find gas masses of 2.4 
xlQi^ 6.9 xlO^ and 1.3 xlO^° M0(each with - 20% 
statistical uncertainty, but note discussion of systematic 
uncertainties above) for Q1623-BX453, Q0449-BX93, 
and DSF2237a-C2 respectively. The stellar mass of 
Q1623-BX453 is 3.1 xlO^" Mg, for a gas fraction /i = 
0.44; for DSF2237a-C2 we find = 1.6 x 10^° Mgand 
/i = 0.32. We are unable to determine a gas fraction 
for Q0449-BX93 because of the lack of near-IR imaging 
from which to derive a stellar mass, although the mass 
limit estimated from the mass-metallicity relation (Erb 
et al. 2006c) suggests that /i > 0.41. Combining the typ- 
ical uncertainties of the stellar (ctm* /M* = 0.4) and gas 
('''Mgaa /Mgas = 0.2) masses, we estimate the uncertainty 
in the gas fraction (JpZ/i — 0.35. The stellar mass and 
gas fraction of Q1623-BX453 are close to the mean val- 
ues (< Ah >= 3.6 ±0.4 X 10i°Mo, < n >= 0.37 ±0.03) 
found for a larger sample of 114 galaxies at z ~ 2 by Erb 
et al. (2006b). While there is no similar determination 
of the mean gas fraction at z ^ 3, the stellar mass of 
DSF2237a-C2 is only slightly higher than the median of 
1.4 X 10^° Mg found for z ~ 3 Lyman break galaxies by 
Shapley et al. (2001; we convert their reported median 
to our adopted IMF). 

We also use the extinction-corrected SFRs determined 
above to estimate the average and peak star formation 
rate surface density Ssfr (Table [7]). Dividing the SFRs 
by the total areas, we find an average Esfr = 8, 4, and 
15 Mq yr-i kpc-2 for Q1623-BX453, Q0449-BX93 and 
DSF2237a-C2 respectively. These values lie in the middle 
of the distribution of Ssfr found for local starbursts by 
Kennicutt 1998 (after accounting for the different IMF), 
and at the upper end of the range for z ~ 2 galaxies 
found by Erb et al. (2006b). Based on the peak flux 
densities observed in a given spaxel, we estimate the cor- 
responding peak star formation rate surface densities to 
be roughly twice the mean, with Esfr — 20, 7, and 25 
Mq yr^^ kpc^^. These values are comparable to the 
maximal star-formation rate densities observed in the lo- 
cal universe by Lehnert & Heckman (1996a; Eiim = 11 
Mq yr^^ kpc^^ in our adopted IMF) and Meurer et al. 
(1997; Eiiin = 25 M© yr^^ kpc^^ in our adopted IMF), 
implying that Q1623-BX453 and DSF2237a-C2 may be 
forming stars at close to the maximal rate in their bright- 
est regions. Given that (in the local universe) a SFR 
density of = 0.1 A/q yr^^ kpc^^ is sufficient for the 
resulting supernovae to drive a wind into the surrounding 
ISM (Heckman 2002), it is likely that rapidly star form- 
ing galaxies similar to these are quite "leaky" systems, 
with a significant fraction of metal-enriched gas finding 
its way into the IGM (Erb et al. 2006c; Adelberger et al. 
2003,2005b). 



6. 



SPATIALLY RESOLVED KINEMATICS AND 
DYNAMICAL MASSES 



Using the integrated spectra described in §4.1 we es- 
timate the kinematic dispersion of the ionized gas by 
fitting a gaussian profile (corrected for instrumental res- 
olution) to the observed emission. We obtain values of 
av = 92 ± 8, 72 ± 12, 101 ± 15 km s'^ respectively for 
Q1623-BX453, Q0449-BX93, and DSF2237a-C2. 



Although our velocity dispersion for Q1623-BX453 is 
slightly higher than that obtained previously by Erb et 
al. (2006b; CTv = 61 ± 4 km s''^)^^, the measured disper- 
sion of DSF2237a-C2 agrees with the value ctv = 100 ± 4 
km s~^ observed by Pettini et al. (2001). However, 
as we demonstrate below in §6.3 the velocity shear of 
DSF2237a-C2 contributes to the line width of the inte- 
grated spectrum. Based on examination of the spatially 
resolved velocity map of the galaxy, we conclude that the 
central velocity dispersion CTv = 79 km s~^, rising slightly 
to about 89 km s^^ at larger radii. 

While the total amplitude of any velocity shear may 
be an unreliable probe of the dynamical mass of a sys- 
tem, the central velocity dispersion a is predominantly 
driven by the gravitational potential of the system even 
in kinematically complex systems such as local ULIRGS 
(Colina et al. 2005). These a may therefore be used to 
estimate the dynamical mass of the galaxies using the 
relation 

Mdyn = (4) 

where a relatively straight-forward derivation (e.g. Erb 
et al. 2006b) gives the constant geometric prefactors 
C ~ 5 appropriate for a uniform sphere and C = 3.4 
for a thin disk with average inclination. Based on the 
kinematic data presented in §6.1 - 6.3, we adopt C — 5 
for Q1623-BX453 and Q0449-BX93, and C = 3.4 for 
DSF2237a-C2. 

While a is reasonably well-determined, r is not and 
can affect the resulting dynamical mass calculated via 
Eqn. |4] significantly. Ideally, we might adopt the virial 
radius, for which the derived Mdyn should be a reason- 
able estimate of the total gravitational mass. However, 
the only observationally-determined length scale is that 
of the much smaller region from which nebular line emis- 
sion is detected (see Table [2]). Adopting this nebular 
emission radius, we calculate dynamical masses of Mdyn 
= 17 ± 2, 7 ± 2, and 5 ± 1 xlO^M© for Q1623-BX453, 
Q0449-BX93, and DSF2237a-C2 (using = 79 km s^^) 
within a radius r — 1.74, 1.23, and 1.03 kpc respectively. 
Unsurprisingly (given that our dynamical masses probe 
only the central 1-2 kpc), these masses differ by ~ a fac- 
tor of 15 from the typical halo mass within the virial 
radius calculated for similar galaxies by Adelberger et 
al. (2005a) based on clustering statistics. A previous es- 
timate of M(jyn has also been obtained for Q1623-BX453 
by Erb et al. (2006b); their value Mdyn = 12 x IO^Mq 
is similar to our own, with the predominant difference 
arising from the difference in observed a^- 

We construct kinematic maps of the three target 
galaxies by fitting gaussian profiles to the spectra con- 
tained within each resampled spaxel. Generally, single- 
component models suffice to represent the emission line 
profiles (although see discussion in §6.2 for Q0449-BX93). 
Such maps of the relative radial velocity and velocity dis- 
persion (Fig. [U center and right panels respectively) are 
intrinsically noisy; we reduce the contribution of noisy 
spaxels to the resulting maps by considering only those 

This discrepancy appears to be due to the presence of a resid- 
ual OH sky feature in the NIRSPEC Ha spectrum; the OSIRIS 
velocity dispersion is roughly in agreement with the NIRSPEC ve- 
locity dispersion of CTv = 106 ±8 km s~^ measured using the [Olll] 
emission line. 
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for which the derived velocity is within 300 km of the 
systemic redshift, the FWHM is greater than or equal to 
the instrumental width, and the SNR of the detection is 
greater than six. These criteria are not derived from 
theoretical limits, but rather represent empirical con- 
straints which we have found (based on iterations with 
various parameters) to produce the cleanest kinematic 
maps. These velocity maps are generally robust to sys- 
tematic uncertainties (since all spaxels are processed in 
a similar manner) and the largest source of uncertainty 
is from random errors arising from the difficulty in de- 
termining the accurate centroid for an emission line of 
given S/N ratio. Given the OSIRIS spectral resolution 
R ~ 3600, the velocity FWHM is about 80 km s'^. We 
estimate that uncertainties in the velocity maps are '^10 
km s^^ for bright central regions of a galaxy, tapering to 
~ 20 km s^^ towards the edges. 

6.1. Q1623-BX453 

The 2d velocity map of Q1623-BX453 (Fig. [J top mid- 
dle panel) is largely consistent with a flat velocity fleld 
across the entire galaxy (~ 4 kpc along the major axis) 
and does not present any evidence for spatially resolved 
shear. Typical deviations from zero relative velocity are 
within the uncertainty of our ability to accurately mea- 
sure velocity centroids from individual spaxels (i.e. about 
10-20 km s~^). The only suggestion of kinematic sub- 
structure is a small ~ 100 mas (800 pc) feature at the 
northeastern edge of the galaxy which appears to be re- 
ceding at ~ 50 km s~^ relative to systemic and has a 
sharp discontinuity with the bulk velocity field (i.e. it 
does not appear to "blend in" by means of a gradual ve- 
locity transition). Given the extremely low flux of this 
feature (which lies outside the faintest flux contour illus- 
trated on Fig. [T]) this may be simply a residual noise 
artifact. 

Although the line-of-sight velocity of the ionized gas is 
roughly constant across the face of the galaxy, the dis- 
persion of this gas (cr) varies slightly with position (Fig. 
m top right panel). In particular, a peaks around 100 
km s~^ in the northeastern region approximately coinci- 
dent with the flux peak, decreases to ^ 70 km s^^ in the 
southwestern region, and drops to below 30 km s~^ in the 
extremely faint emission region to the west. While there 
appear to be peaks of cr ~ 120 km on the northern 
and eastern edges of the galaxy, it is uncertain whether 
these represent real features or are simply noise in the 
extremely faint surface brightness edges of the galaxy. 

The flat velocity field of this galaxy is particularly in- 
triguing in light of the semi-analytic picture of galaxy 
formation in the early universe (e.g. Mo, Mau, & White 
1998, Baugh 2006) in which (in brief) cold gas falls into 
collapsed dark halos, shock heats, radiatively cools to 
form a disk supported by residual angular momentum, 
and proceeds to form stars within this gaseous disk. Un- 
less Q1623-BX453 happens to be inclined almost per- 
fectly face-on to our line-of-sight^^, the absence of a ve- 
locity gradient across the galaxy indicates that the kine- 
matics of the ionized gas are not dominated by ordered 

^® In order to mask a 100 km circular velocity within a 20 
km uncertainty a galaxy must have inclination i < 10°. This 
occurs in less than 2% of cases for a randomly distributed set of 
inclinations. 



rotation. Clearly, the simplest description of galaxy for- 
mation is incomplete. 

One variant on the basic rotating disk model which 
provides a reasonable explanation of some of the ob- 
served characteristics is that described by Noguchi et al. 
(1999) and Immeli et al. (2004; their models A/B), in 
which highly-efficient cooling mechanisms lead to early 
fragmentation of the gas disk into self-gravitating clouds. 
These clouds are sufficiently massive that they spiral to 
the center of the galaxy via dynamical friction within a 
few dynamical times, largely dragging the existing stellar 
population along with the gas. Immeli et al. (2004) find 
that an intense burst of star formation typically occurs in 
the nuclear region of the galaxy when these massive gas 
clouds collide, typically peaking at around 100 Mq yr~^ 
(over twice that expected for "traditional" star formation 
histories in galaxies with less efficient gas cooling mech- 
anisms) and leading to the rapid formation of a stellar 
bulge. Such a theory may provide a natural explanation 
for Q1623-BX453, which we infer to have a large cold 
gas fraction, near-peak star formation rate density, and 
a sizeable stellar population. Alternatively, such char- 
acteristics may also be explained by the merger of two 
massive, gas-rich galaxies whose cold gas reservoirs are 
undergoing rapid star formation in the center of the new 
gravitational potential. 

While it might appear surprising that a system with 
such violent dynamics should show so little evidence for 
kinematic structure, such systems are well-known in the 
local universe. Local ULIRGS may often have gas masses 
constituting the bulk of their dynamical mass (Solomon 
et al. 1997), fueling intense starbursts from the dense 
molecular gas in their cores. The common presence of 
multiple nuclei and tidal features (e.g. Sanders et al. 
1988, Bushouse et al. 2002) suggests that these star- 
bursts may commonly arise from the major merger of two 
gas-rich systems, and the ionized gas kinematics typically 
show star formation in gas flows rapidly forming stars 
in the center of the gravitational potential. The veloc- 
ity fields of ULIRGS are typically dominated by strong, 
asymmetric gas flows (Colina et al. 2005), although the 
specific properties of these flows vary considerably from 
case to case. The velocity structure of Q1623-BX453 
resembles in particular the kinematics of Mrk 273 and 
IRAS 15250-f3609 (see Fig. 1 of Colina et al. 2005); 
these galaxies show no strong kinematic substructure in 
the brightest regions of Ha emission within ^2 — 3 kpc 
of the peak Ha emission despite relatively strong (~ 200 
km s~^) kinematic features in regions with Ha surface 
brightness lower by a factor ~ 5. It is an interesting spec- 
ulation therefore whether deeper observations of Q1623- 
BX453 will reveal similar kinematic structure in regions 
of lower surface brightness. 

6.2. Q0U9-BX93 

While the integrated spectrum of Q0449-BX93 (Fig. 
[2]) is dominated by a single emission spike deflning the 
systemic redshift, the spectrum of the northwestern re- 
gion of the galaxy shows a faint secondary emission peak 
blueshifted by ~ 180 km s""'^ from the systemic red- 
shift (Fig. [6|). Given that the FWHM of the primary 
emission component is 150 km s~^, these two com- 
ponents in Figure [S] are blended, motivating the use of 
a 2-component gaussian model in order to adequately 
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characterize the emission hne profile of the source. 

Considering first the primary spectroscopic compo- 
nent, the velocity and dispersion structure of Q0449- 
BX93 (Fig. [U center and middle-right panels respec- 
tively) are similar to those of Q1623-BX453 although 
there is mild evidence for velocity shear along the 
northwest-southeast axis. The peak-to-peak amplitude 
of this velocity differential is only ~ 40 km s~^, and 
therefore not highly significant with respect to our es- 
timated velocity uncertainty of 10 - 20 km s~^. Even 
if this shear is genuine, it is not obvious that this im- 
plies a rotationally supported system. Assuming a sim- 
plistic model in which the circular velocity at the outer 
edge of the galaxy (1.23 kpc; see Table [J) is 40 km 
(i.e. with an extremely generous inclination correction 
to the observed 20 km s~^), we obtain a dynamical mass 
of only Mdyn = 5 x lO^M©, far less than the gas mass 
alone within the same region as derived from the Ha 
emission flux. That is, the observed shear is insufficient 
to provide dynamical support for the galaxy. As in the 
case of Q1623-BX453, the velocity dispersion is consid- 
erably larger {a — 72 km s^^) than the peak-to-peak ve- 
locity shear, suggesting the physical interpretation that 
Q0449-BX93 may also represent a merger-induced star- 
burst or the bright central emission region of a gas-rich 
disk which has fragmented due to instability. While in- 
sufficient broadband data exist for Q0449-BX93 to be 
able to reliably fit stellar population models, estimates 
of the stellar mass using the mass-metallicity relation of 
Erb et al. (2006c) suggest that this galaxy is consider- 
ably less massive (Af* < 10^°Mq) than Q1623-BX453. 
Coupled with its moderate extinction-corrected Ha star 
formation rate of 18 Mq yr~^ (close to the average for 
the z - 2 galaxy sample, Erb et al. 2006a) Q0449-BX93 
may be more representative of the general galaxy popu- 
lation. 

Selecting next the gaussian fits to the secondary emis- 
sion line component it is possible to create velocity and 
dispersion maps (Fig. [7]) of this faint, kinematically dis- 
tinct sub-component (hereafter KDSC). To the limiting 
surface brightness of our study, the KDSC is marginally 
resolved over only a small (~ 200 x 150 mas, i.e. ~ 1 
kpc) part of the northwestern region of the galaxy com- 
parable in size to the PSF (~ 150 mas). The kinematic 
offset of the KDSC from the systemic redshift is approx- 
imately constant across the region of emission and does 
not smoothly join to the kinematics of the primary com- 
ponent. While this may be simply due to the small size 
of the KDSC (across which we expect only ^1 — 2 spa- 
tially independent spectra), an additional interpretation 
(within the framework of a fragmenting disk model) may 
be that the KDSC represents one of multiple starburst- 
ing regions in the process of coalescing in the center of 
the galaxy for which we happen to see a kinematic dif- 
ferential. Alternatively, the KDSC may have initially 
begun star formation as a unit (i.e. satellite galaxy) dis- 
tinct from the bright component of Q0449-BX93 which 
is only now merging with its host galaxy and is perhaps 
responsible for the observed burst of star formation. It 
is unfortunately beyond the scope of the current data 
to derive the parameters (e.g. relative mass) of such a 
scenario with any pretense of reliability. 

6.3. DSF2237a-C2 




Fig. 6. — The spectrum of Q0449-BX93 is shown (in units of 
velocity relative to systemic) integrated over the region in which 
a secondary Ha emission component is observed blueward of the 
systemic emission. Vertical dashed lines denote the fiducial galaxy 
velocity (i.e. km s~^) and the location of the secondary emission 
peak after deblending. This peak occurs at a relative velocity of 
-177 ±11 km s-i. 



The kinematic map of ionized gas in DSF2237a-C2 is 
notably different from those of Q1623-BX453 and Q0449- 
BX93, presenting clear evidence for resolved velocity 
shear along the morphological major axis at redshift 
z > 3. As illustrated in Figure [T] (bottom center panel), 
the velocity of peak [O ill] emission relative to the sys- 
temic redshift varies smoothly from about —70 to +70 
km s^^ over a distance ~ 2.5 kpc. An apparent ridge-line 
of high velocity (most evident in the southwest corner of 
the galaxy) presents some evidence for a fiattening in the 
velocity curve, although this is far from robust. The ve- 
locity dispersion of individual regions within the galaxy 
is (7 ~ 79 km s~^; this is inflated by velocity shear in the 
integrated spectrum, giving a — 101 km s^^. 

For comparison to previous studies we produce an 
"ideal" long-slit spectrum of DSF2237a-C2 by extract- 
ing spaxels along the major axis (PA 38°, see §3) us- 
ing a 180 mas wide "slit" and plot the resulting velocity 
curve in Figure [51 The gradual trend of velocity across 
the galaxy is clear in this figure, along with the appar- 
ent fiattening of the velocity curve in the outskirts of 
the galaxy. Although long-slit spectra were previously 
obtained with Keck-NIRSPEC observations (Pettini et 
al. 2001), no shear was observed by these authors de- 
spite close alignment of the NIRSPEC slit (PA 55°) with 
the now- known shear axis (PA 38°). Given that the to- 
tal peak-to-peak amplitude of the shear (~ 140 km s~^) 
should be identifiable with the spectral resolution of NIR- 
SPEC, we infer that DSF2237a-C2 represents a practi- 
cal case-example of how seeing-limited studies may miss 
velocity structure on spatial scales smaller than that of 
the seeing halo. For comparison, in Figure [9] we plot 
the 2-d velocity map recovered from our OSIRIS obser- 
vations of DSF2237a-C2 after introducing a 0.5" spatial 
gaussian blur to roughly simulate the smearing effects 
of atmospheric seeing. While the overall semblance of 
shear is retained, the peak-to-peak observed amplitude 
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Fig. 7. — As Figure[T] shows the Ha emission hne velocity and velocity dispersion maps for the faint blucshiftcd component of Q0449-BX93. 
Overlaid contours are Ha flux density. 



Distance (kpc) 

Fig. 8. — One-dimensional spectrum of DSF2237a-C2 (solid 
boxes) extracted from the OSIRIS data using an idealized "slit". 
Uncertainties are conservative estimates based on the number of 
spaxels (each with 20 km s~^ uncertainty) included in the spec- 
trum at each position along the simulated slit. The solid line indi- 
cates the FWHM of the PSF. Note that this figure underestimates 
the turnover velocity apparent in 2-d data because spatial aver- 
aging across the simulated slit combines peak ridgeline velocities 
with lower velocities on either side of the "v" -shaped pattern (see 
Fig. [lOj. 

of this shear is reduced by a factor ~ 2.^^ In concert 
with the lower spectral resolution of NIRSPEC (in the 
low-resolution mode used for the Pettini et al. 2001 ob- 
servations) and slight offsets in slit placement and orien- 
tation from the optimal, it is easy to see why previous 
observations did not detect shear. 

Adopting the inclination and position angle derived in 
§3 from the nebular morphology, we create an idealized 
thin-disk model (i.e. a rotating sheet) which reproduces 
the size and orientation of DSF2237a-C2 and adjust the 
asymptotic velocity and turnover radius as needed to 
match the observations shown in Figures [T] and [8l This 

^'^ We note that Erb et al. (2004; their Fig. 3) reach a simi- 
lar general conclusion based on observations of the galaxy Q1700- 
BX691, for which the apparent amplitude of velocity shear varied 
by at least a factor of ~ 2 depending on the quality of seeing. 




Fig. 9. — Velocity map recovered from OSIRIS observations of 
DSF2237a-C2 after convolution with a 0.5" PSF. The field of view 
is identical to that in Figure [T] 

model is characterized by a radius of 0.21" = 1.6 kpc and 
a deprojected asymptotic velocity of 93 km s~^, implying 
a dynamical mass of 3.2 x lO^M© interior to 1.6 kpc. 
We plot the velocity map of this model disk in Figure [TO] 
(left panel) , along with the velocity map recovered from 
this model using simulated OSIRIS observations which 
include convolution with the PSF and extraction of the 
high surface brightness pixels (right panel). We have not 
attempted to match every detail of the observed velocity 
map by fitting a more complex model, and some discrep- 
ancies remain (the blue-side velocities in particular are 
imperfectly matched by the simple model). While more 
complicated fitting routines could certainly be performed 
(and indeed have been for galaxies studied by Forster 
Schreiber et al. 2006 and other groups), the exacting fits 
of such models frequently belie the large uncertainties in- 
herent to the observations. For example, the evidence for 
a turnover in the velocity curve of DSF2237a-C2 largely 
disappears with only slightly less careful data reduction, 
and the salient features of the model are almost indistin- 
guishable from similar models for which the inclination 
of the disk to the line of sight differs by over 20°. 

Given the uncertainties in such quantitative models, 

This mass is slightly lower than that calculated on the basis 
of the velocity dispersion (and quoted in Table [Sjl, refiecting in 
large part the uncertainty in the geometric prefactors used when 
converting from velocities to estimates of the dynamical mass. 
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what may we conclude about this galaxy with some con- 
fidence? It is unlikely that the shear in DSF2237a-C2 
is due primarily to a major merger because the veloc- 
ity field shows a smooth, gradual transition rather than 
any sharp discontinuity as might be expected for two 
kinematically distinct components. In addition, the flux 
distribution is smooth and centrally concentrated; we do 
not detect multiple nuclei on scales more widely sepa- 
rated than ^ 750 pc (the spatial resolution of the obser- 
vation). Adopting the simplest interpretation that the 
observed shear represents a rotating gaseous disk, it is 
worth asking whether the disk is stable against gravita- 
tional collapse. One traditional way of quantifying such 
stability is via the Toomre 'Q' parameter (Toomrc 1964) 
which effectively assesses the mass of the disk with re- 
spect to the total gravitating mass of the system: 

LrAidisk/'Tdisk JWdisk 

Typically, values of Q < 1 represent dynamically unsta- 
ble systems for which the total baryonic mass exceeds 
that which can be supported by the observed rotational 
velocity of the disk, indicating either that the disk itself 
is unstable or that non-rotational motions contribute a 
significant degree of support to the system. In contrast, 
values Q > 1 suggest that the observed baryonic compo- 
nent can be supported by the observed rotational motion 
alone, while Q ^ 1 represents the case of marginal sta- 
bility. 

Considering all of the uncertainties involved in deter- 
mining the baryonic and dynamical masses, it is not pos- 
sible to establish a single well-defined value for Q, al- 
though realistic estimates for various cases place a limit 
Q < 1. In the first case, we consider the ratio of masses 
within the radius probed by the OSIRIS detections (r = 
1.03 kpc), for which Mdisk = + = 4.2 x IO^^Mq 
in the limit that the entirety of the stellar mass lies within 
this radius, and Mdisk = 1.3 x lO^^M© in the opposite 
limit for which there is negligible stellar mass within the 
central kpc. Similarly, we may elect to use estimates of 
the dynamical mass based on the central velocity disper- 
sion (Mdyn = 5 X 10^ Mq) or the rotationally-supported 
thin disk model (Mdyn = 3.2 x lO^Mg). Combining 
these assorted estimates, we obtain a range of values 
0.08 < Q < 0.38. In the second case, we calculate the 
global Q value if we consider the mass within a much 
larger 10 kpc radius. At such radii Mdisk > 4.2 x lO^^M© 
(as the contribution of gas mass outside the central kilo- 
parsec is unknown). The dynamical mass is even more 
uncertain^^: inserting r = 10 kpc in Equation |4] we find 
-^dyn = 5 X lO^^M©, while extending the disk model 
with asymptotic velocity 93 km s^^ out to such radii 
gives Mdyn = 2 X 10^'^ Mq. Adopting the highest physi- 
cally motivated limit, we determine that Q < 1.2 within 
10 kpc. In both cases (and most noticeably within the 
central kpc), the baryonic mass in the system outweighs 
what could be supported by a simple rotating disk model, 
highlighting the importance of the ionized gas velocity 
dispersion and indicating the improbability that star for- 
mation in DSF2237a-C2 occurs predominantly in a rota- 

The dynamical mass may be particularly uncertain if the ve- 
locity curve continues rising out to larger radii instead of turning 
over at ^ 93 km s^^ as suggested by observations. 



tionally supported thin disk of gas (see also discussion by 
Forster Schreiber et al. 2006). Interestingly, the kine- 
matics of this galaxy are quite similar to those of the 
z = 1.4781 galaxy Q2343-BM133 observed by Wright et 
al. (2007), which also had dispersion (cr = 73 ± 9 km 
s^^) comparable to its velocity shear (one-sided ampli- 
tude V — 67 km s~^). 

7. DISCUSSION 

We have presented the results of a pilot study of the 
spatially resolved spectra of galaxies at redshift z ~ 2 — 3, 
demonstrating that the velocity structure of the ionized 
gas appears to be inconsistent with simple rotational 
support mechanisms. While there has been some ef- 
fort to target galaxies with particularly high star for- 
mation rates to ensure detection in these early stages 
of the observational program, the Ha flux of one of our 
three galaxies (Q0449-BX93, whose star formation rate 
was predicted from UV photometry) is found to be quite 
close to the average (Erb et al. 2006a, b) of the z ^ 2 
optically-selected galaxy sample. 

It is, nevertheless, more than a little challenging to 
attempt to draw general conclusions about the entire 
z 2 — 3 galaxy population on the basis of only three 
galaxies, and therefore interesting to compare our re- 
sults to integral-field velocity maps obtained recently by 
Forster Schreiber et al. (2006) for 14 rest-UV selected 
galaxies at redshift z ^ 2 drawn from the catalog of 
Erb et al. (2006a) and observed with SINFONI (Eiscn- 
hauer et al. 2003) on the VLT. While selection criteria 
varied, generally galaxies were known to be large and 
bright and (in many cases) known to exhibit consider- 
able velocity shear or dispersion on the basis of previous 
long-slit observations (Erb et al. 2003, 2006a). Of these 
14 galaxies, 9 show evidence of spatially resolved veloc- 
ity gradients with peak-to-peak amplitude ranging from 
40 to 410 km s~^ over ~ 10 kpc. Certainly, systems 
with strong kinematic features appear to exist at z ~ 2, 
some of which may genuinely represent orbital veloci- 
ties in massive rotating disks. Indeed, the mean stellar 
mass of the subsample with shear greater than 50 km 
s^^ is 9.4 X IO^^Mq (based on the 6 galaxies for which 
M^ is given by Erb et al. 2006b), almost three times 
greater than the mean of the Erb et al. (2006b) sample 
(< >= 3.6 ±0.4 X IO^^Mq). In contrast, the remain- 
ing galaxies in the Forster Schreiber et al. (2006) sample 
which do not show such large-scale shear have mean stel- 
lar mass 2.6 x lO^^M©, similar to the masses of Q1623- 
BX453, Q0449-BX93, and DSF2237a-C2 (M* = 3.1, < 1, 
and 2.9 x lO^^M© respectively). Rather than represent- 
ing rare cases (such as counter-rotating mergers; Forster 
Schreiber et al. 2006), the distribution of stellar masses 
suggests that these non-rotating galaxies may instead be 
quite common in the z ~ 2 population. 

Unfortunately, it is difficult to directly compare the de- 
tailed kinematics of the Forster Schreiber et al. (2006) 
observations with our own since these data are seeing- 
hmited with PSF FWHM ~ 0.6" (therefore smearing 
kpc-scale kinematics like those discussed in this contri- 
bution beyond recognizability; see Fig. [5]). A more de- 
tailed comparison may be made to the z = 2.3834 galaxy 
BzK-15504 observed by Gcnzel et al. (2006) using SIN- 
FONI in combination with natural guide-star adaptive 
optics for an angular resolution of 150 mas (i.e. com- 
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SPECTROSCOPY OF HIGH REDSHIFT GALAXIES 

DSF2237a'^C2 Madel Disk 




Fig. 10. — Left panel: Idealized rotating disk model for DSF2237a-C2. Right panel: Velocity field which might be observed with OSIRIS 
after convolution of the idealized disk model with the real PSF and selection of an appropriate distribution of bright "emission" regions. 
Compare to the actual recovered velocity map shown in Figure [T] (bottom center panel). The field of view is identical to that in Figure [T] 



parable to the ~ 110 - 150 mas resolution of our target 
galaxies after smoothing). While this galaxy exhibits re- 
solved velocity shear in the outer regions of the galaxy 
it is not well fit by a simple disk model within 3 kpc 
of the nucleus (i.e. approximately the region probed in 
our three target galaxies), which Genzel et al. (2006) 
postulate may be due to radial gas flows fueling the cen- 
tral AGN whose presence is inferred from rest-frame UV 
spectroscopy. One possibility may therefore be that dis- 
ordered nuclear kinematics could generally join to more 
orderly shear at larger radii, although in our data we find 
evidence for neither AGN nor large-radius shear. Indeed, 
even in the circum-nuclear regions of BzK-15504 it is not 
obvious that a simple disk model provides an optimal de- 
scription of the kinematics since the line-of-sight velocity 
dispersion remains high at large radii (cr ~ 60 — 100 km 
s~^, similar to the large values observed in DSF2237a- 
C2, and in Q2343-BM133 by Wright et al. 2007), and 
difference maps of the observed velocity field with the 
best-fitting exponential disk model show deviations of 
greater than 100 km s^^ (see discussion by Wright et al. 
2007). 

What, then, may we conclude about the kinematics 
of these early galaxies and their implications for galaxy 
formation? One common explanation for both the un- 
usual kinematics and the high nuclear gas fraction is that 
these galaxies may represent major mergers of gas-rich 
galaxies. While this explanation may suffice for a few 
cases however (e.g. perhaps Q0449-BX93), it is unlikely 
(based on current cosmological models) that this is the 
general explanation for the z ^ 2 population since the 
space density of such galaxies with similar physical prop- 
erties is far too high. Alternatively, efficient gas cool- 
ing in the early universe could have led to high cold gas 
fractions in early-stage galactic disks which then frag- 
mented due to self-gravity and collapsed to form a nu- 
clear starburst (e.g. Immeli et al. 2004a). In both pic- 
tures, considerable kinematic structure should exist in 
these galaxies although (similar to local ULIRGS; Col- 



ina et al. 2005) the kinematics of the brightest nebular 
emission regions may be relatively featureless and quite 
variable in detail from system to system. However, it is 
unknown whether the lack of shear (or indeed any signifi- 
cant velocity structure beyond a high dispersion) persists 
to lower surface-brightness ionized gas at larger radii or if 
the kinematics of this gas are strong and disordered, or- 
dered with dynamically cold orbital motion, or if indeed 
significant nebular line emission exists at all beyond the 
radii probed in this study. For one of the galaxies studied 
here (Q1623-BX453) previous deep long-slit spectroscopy 
(Erb et al. 2006b) could have detected any significant 
shear if it were present at larger radii (and the kinematic 
major axis was reasonably well-aligned with the slit), but 
found no evidence for such a trend. 

What is perhaps most clear is that we do not yet under- 
stand the dynamical state of galaxies during this period 
in which they are forming the bulk of their stars. While 
traditional theories of galaxy formation posit that the 
majority of star formation should occur in rotationally 
supported gaseous disks, this does not appear to be in 
good agreement with the accumulating observational ev- 
idence, inviting the construction of new models of gas 
cooling which can explain the flat, dispersion-dominated 
kinematic features of galaxies in the young universe. 
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